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ABSTRACT 

We present the determination of the geometric i?-band albedos of two main-belt 
comet nuclei based on data from the Spitzer Space Telescope and a number of ground- 
based optical facilities. For 133P/Elst-Pizarro, we find an albedo of pr = 0.05 ± 0.02 
and an effective radius of r e = 1.9 ± 0.3 km (estimated semi-axes of a ~ 2.3 km and 
b ~ 1.6 km). For 176P/LINEAR, we find an albedo of p R = 0.06 ±0.02 and an effective 
radius of r e = 2.0±0.2 km (estimated semi-axes of a ~ 2.6 km and b ~ 1.5 km). In terms 
of albedo, 133P and 176P are similar to each other and are typical of other Themis family 
asteroids, C-class asteroids, and other comet nuclei. We find no indication that 133P 
and 176P are compositionally unique among other dynamically-similar (but inactive) 
members of the Themis family, in agreement with previous assertions that the two 
objects most likely formed in-situ. We also note that low albedo (pr < 0.075) remains 
a consistent feature of all cometary (i.e., icy) bodies, whether they originate in the inner 
solar system (the main-belt comets) or in the outer solar system (all other comets). 

Subject headings: comets: general — minor planets, asteroids 



lr This work makes use of observations made with the Spitzer Space Telescope (Programs 3119 and 30678), which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with the National 
Aeronautics and Space Administration (NASA). Additionally, some data presented herein were obtained at the W. 
M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, 
the University of California, and NASA, and was made possible by the generous financial support of the W. M. 
Keck Foundation. Some data presented herein were also obtained at ESO facilities at La Silla under program ID 
081.C-0822(A). 
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INTRODUCTION 



The main-belt comets (MBCs), of which 133P/Elst-Pizarro (hereafter, 133P) and 176P /LINEAR 
(hereafter, 176P) are examples, occupy stable orbits that are de coupled from Jupiter and which 
are indistinguishable from the orbits of other main-belt asteroids (IHsieh &: Jewittll2006bl ). Dynam- 
ical simulations show that MBCs are extremely unlikely to originate in the Kuiper Belt given the 
current configuration of the major planets (e.g.. iFernandez et a/.ll2002j ). indicating that they are 
instead likely to be native to the main asteroid belt. Recent work suggests that some icy Kuiper 
Belt objects might have been delivered to the asteroid belt during the Late Heavy Bombardment 
(jLevison et aZ.ll2Q08l ) . but even those simulations fail to produce the low-inclination, low-eccentricity 
orbits of MBCs such as 133P and 176P. 



In this letter, we use observations from the Spitzer Space Telescope (hereafter, Spitzer; IWerner et al. 



2004) to determine the geometric albedos of 133P and 176P, and then discuss the implications of 
these measurements. 



2. OBSERVATIONS 



We obtained optical observations of 133P and 176P on multiple occasions from 2003 through 
2008 using the 10 m Keck I and University of Hawaii (UH) 2.2 m telescopes on Mauna Kea, and 
the 3.58 m New Technology Telescope (NTT) at the European Southern Observatory (ESO) at La 
Silla. Observations with the UH 2.2 m telescope were made using either a Tektronix 2048 x 2048 
pixel CCD or the Orthogonal Parallel Transfer Imaging Camera (OPTIC; iTonrv et aZ.ll2004l ) . both 
behind standard Kron-Cousins BVRI broadband filters. Observa tions with Keck were made using 
the Low Resolution Imaging Spectrometer (LRIS: lQke et aZ.lll995l ) in imaging mode. LRIS employs 
a Tektronix 2048 x 2048 CCD with standard Kron-Cousins BVRI filters. Observations with the 



NTT were made using the ESO Faint Object Spectrograph and Camera (EFOSC2; iBuzzoni et al 



19841 ). which employs a 2048x2048 pixel Loral/Lesser CCD behind Bessel BVR broadband filters. 



Bias subtraction and flat-field reduction were performed for all optical data. Dithered images of 
the twilight sky were used to construct flat fields for UH 2.2 m data, while images of the illuminated 
interior of the telescope dome were used to construct flat fields for Keck and NTT data. Photometry 
of our target objects and lLandoltl (|1992l ) standard stars was obtained by measuring net fluxes within 
circular apertures of varying radii depending on the nightly seeing, with background sampled from 
surrounding circular annuli. 

Spitzer observations of 133P (3 visits, 166 s of total exposure time per visit; Fig. [lh.), using the 
24 /im channel (effective wavelength of 23.68 |im) of the Multiband Imaging Photometer for Spitzer 
(MIPS:lRieke et a/Jl2004r ) and originally obtained on 2005 April 11 as part of Cycle 1 program 3119 



(Reach et al 



20071 ). were retrieved from the Spitzer archive. Observations of 176P (2 visits, 48 s 



of total exposure time per visit; Fig. QJ)), also with the 24 fim channel of MIPS, were obtained 
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on 2007 January 1 as part of Cycle 3 program 30678. Observational circumstances are shown in 
Table [TJ Photometry of our target objects from pipeline-processed Spitzer post-Basic Calibrated 
Data (PBCD) was obtained by measuring net fluxes within circular apertures with 6-pixel (14'.'7) 
radii, and then applying appropriate aperture corrections (1.14 in the case of a 6-pixel aperture) 
and color corrections (0.96 for both targets). 



3. RESULTS 

We use our optical data to find best-fit IAU phase function parameters for 133P of Hr = 
15.49 ± 0.05 mag and G = 0.04 ± 0.05, and best-fit linear phase function parameters (omitting 
data obtained at solar phase angles at which opposition surge effects are expected) of m^(l, 1, 0) = 
15.69 ± 0.05 mag and (5 = 0.049 ± 0.004 mag deg -1 . These parameters are calculated using 
photometry obtained while 133P was observe d to b e inactive, and as such, are a refinement of 



parameters previously derived by iHsieh et all (|2004h from photometry obtained while 133P was 



visibly active. For 176P, we find corresponding parameters of Hr = 15.10 ± 0.05 mag, G = 
0.26±0.05, J71r(1, 1,0) = 15.27±0.05 mag, and f3 = 0.034±0.005 mag deg -1 . These parameters were 
likewise calculated only using photometry obtained while the comet was observed to be inactive. 
Plots of phase function solutions for both objects are shown in Figured From their phase functions, 
we estimate our targets' expected mean optical magnitudes as viewed from Spitzer at the time of 
their observations to be mji = 21.63 mag for 133P and ttir = 20.39 mag for 176P. 

Both objects exhibit significant rotational brightness variations, however, which represent sig- 
nificant sources of uncertainty in the interpretation of our infrared data. A rotation period of 
Pmt = 3.471 hr a nd a lightcurve range of Am = 0.4 mag have been previously found for 133P 



(|Hsieh et aZ.I I20Q4J ) . The rotational properties of 176P are currently poorly constrained. On 2007 
March 21, however, we observed a photometric range for the object of Am ~ 0.6 mag over ~ 4.5 hr, 
suggestin g a rotation per i od of P rn t > 18 hr (assuming a double-peaked lightcurve). This is consis- 



tent with iLicandro et all (|2007al ) who found P ro t > 22 hr 



Fortunately, constraints on the rotational phase of each object can be derived from the infrared 
data. Our second flux density measurement for 176P was 1.4 times larger than the first, implying an 
equivalent increase in visible cross-sectional area, corresponding to a change in visual magnitude of 
Am = —0.37 mag, a significant fraction of the object's inferred optical photometric range. Such a 
large magnitude change indicates that the object was necessarily first observed near the minimum 
and then near the maximum of its lightcurve (consistent with the 4.55 hr interval between the 
two observations). Using these constraints, we are able to adjust our optical brightness estimates 
(Table[T|) accordingly, thereby reducing the effects of rotational phase uncertainty. The three Spitzer 
observations of 133P span only 7 min meaning that the rotational phase is less tightly constrained. 
Given 133P's short rotation period, however, a small decrease in the scattering cross-section is 
still detectable between the first and last Spitzer observation, corresponding to a change in visual 
magnitude of Am = 0.14 mag and leading to the revised optical brightness estimates for 133P in 
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Tabled 



We use the [Harris] (|1998l ) Near-Earth Asteroid Thermal Model (NEATM) to iteratively solve 
for the effective radius, r e , and geometric -R-band albedo, pr, of each object. As with any model, 
NEATM requires a number of assumptions, which, in turn, introduce uncertainties. One such source 
of uncertainty is the phase effect for thermal emission. NEATM treats the effect geometrically, 
calculating it based on the fraction of the Earth-facing hemisphere that is illuminated by the Sun 
at the time of observation. While this effect has been poorly measured and thus poorly constrained 
for large phase angles, the infrared phase coefficient of 0.01 m ag deg -1 that we use here is g enerally 
considered to be appropriate for phase angles a < 30° (c/. iMorrisonl 1 19771 ; lHarrisI Il998l ). Thus, 
given the small phase angles (14° < a < 17°) at which the Spitzer observations were obtained, this 
effect should introduce minimal systematic uncertainty into our calculations. 

A more significant issue is that of the beaming parameter, rj. We lack the minimum number of 
data points needed to constrain rj for either 133 P or 176P, forcing us to assume its value. A Spitzer 
survey of ~ 50 Jupiter-family comet nuclei by iFernandez et all (|2008l ). however, found values of 
0.6 < rj < 1.2, and all were consistent with rj ~ 0.94 ± 0.20. Given the results of this survey 
and assuming that 133P and 176P have low therm al inertias (similar to other comet nuclei, e.g., 



9P/Tempel 1, which has / < 50 W K" 1 m -2 s 1 ' 2 ; iGroussin et all 120071 ). we adopt rj = 1.0 as a 
reasonable assumption for solving for r e and pr. To account for uncertainties in rj, we also perform 
parallel calculations for rj = 0.8 and r/ = 1.2 (Table [2]). 

Thus, assuming an emissivity of e = 0.9, we find r e = 1.9 ± 0.3 km and pr = 0.05 =L 0.02 for 
133P, and r e = 2.0 ± 0.2 km and pr = 0.06 ± 0.02 for 176P. Estimated errors for both objects 
are mainly due to uncertainties in both r\ and rotational phase. Given the observed photometric 
ranges (Am^p ~ 0.40 mag; Ami76p ~ 0.60 mag) and corresponding inferred minimum axis ratios 
([a/6]i33P ~ 1-45; [a/6]i76p ~ 1.74) for each object, we find a ~ 2.3 km and b ~ 1.6 km for 133P, 
and a ~ 2.6 km and b ~ 1.5 km for 176P as our best estimates of the semiaxes of each object. 



DISCUSSION 



We plot histograms showing the albedo (pv) distributions of several solar system body pop- 
ulations of interest in Figure [3l Assuming that both 133P and 176P are approximately spectrally 
neutral (i.e., py « Pr)-, based on 133P's spectral classification as a C- or B-type asteroid and 
176P's classification as a B-type asteroid (jLicandro et al\ l2007al ) . we find that their albedos are 
typical of C-class asteroids (Fig. [3k,) and are also well within the distribution of albedos measured 
for members of the Th emis asteroid family ( Fig. Eb), with which 133P and 176P appear to be 
dynamically associated (IHsieh fc Jewittll2006bl ). 



The Themis family is dominated by C-class asteroids (cf. iFlorczak et aZ.1 11999) . of which a 
substantial fraction (17 of the 39 currently classified members of the family, or ~ 44%) belong to 
the subclass of B-type asteroids. For comparison, B-type asteroids comprise only ~ 15% of the 
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general C-class pop ulation and ~ 5% of all currently classified asteroids (jTholen &: Baruccilll989l ; 
Lazzaro et al\ 12004 ) . In terms of albedo, we find 133P and 176P to be consistent with both C-type 
as teroids and B-type aste roids (cf. Figs. [3^ and[3]i), in agreement with their spectral classifications 
bv lLicandro et al\ (|2007al ). Thus, in terms of both albedo and spectral type, 133P and 176P appear 
to be typical Themis asteroids, supp orting previous speculation that the family might be home to 
more MBCs dHsieh fc JewittlhoOBah . 



Other objects like 133P and 176P that also have orbits considered to be dynamically asteroidal, 
yet have been associated with observed or inferred cometary activity, are also classified as C-class 
objects. One such object is the c ross-listed comet-asteroid 107P/(4015) Wilson-Harrington (C- 
or F-type; iTholen &: Bar ucci 198 9|). Other examples include Geminid meteor stream parent 3200 



Phaethon (B- or F-type; ITholen Barucci 



1989 



Licandro et al 
200d: 



2007 bJ) and its possible fragment, 
Another likely 



Kinoshita et fli.ll2007l ). 



155140 (2005 UD) (B- or F-type; I.Tewitt fc Hsiehl 
fragme nt of 3200 Phaethon, 1999 YC, appears spectrally neutral and is classified as a C-type 
object (jKasuga Jewittl 120081 ) . 



Being below the upper bound of "comet-like" albedos (pr = 0.075) employed bv lFernandez et al 



(|2005l ). the albedos of 133P and 176P are also consistent with those of the nuclei of other ac- 
tive comets (Fig. (3]1). Spectroscopically, comet nuclei exhibit a broad range of colors, with 
both D-type-like and C- t ype-like spec tral reflectivity gradients being found for various comets (cf. 



Fitzsimmons et al. 



1994 



Jewittl [2002j) , and in terms of albedos, the two MBC nuclei we consider 



here are consistent with both spectral types (Fig. [3^ &[3fe). Thus, despite their strong dynamical 
association with main-belt asteroids, we find that 133P and 176 P have surfac es that may be com- 
positionally comparable to other comets. This is consistent with I Jewittl ((2002) who suggested that 
the surface properties of short-period comet nuclei were likely largely due to sublimation-driven 
evolutionary effects and were not primordial in nature. 

Studying the surface properties of the main-belt comets is vital for understanding their evo- 
lution and putting their volatile content into the proper context. In light of those goals, we find 
that, in terms of albedos, (1) 133P and 176P are similar to each other, (2) they are typical of other 
Themis asteroids and the C- and B-type asteroids that dominate the Themis family, and (3) their 
albedos are also consistent with albedos measured for other comet nuclei and D-type asteroids. 
Given these results, we find that low albedo continues to be a consistent feature of all cometary 
bodies, whether they originate in the outer or inner solar system. This finding necessarily also 
means, however, that albedo does not appear to be an effective diagnostic of the region from which 
a comet originates. 



We acknowledge support of this work through STFC fellowship grant ST/F011016/1 to HHH, 
NASA Spitzer grant JPL-1289078 and NASA Planetary Astronomy grant NNG05GF76G to DJ, 
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Table 1. Spitzer Observations 



Object 


Date 




UT 


R a 


A Sp h 




m R 


d 


133P/Elst-Pizarro 


2005 


Apr 


11 


08 


01:11 


3.596 


3.046 


14.6 


21.56 


± 


0.15 




2005 


Apr 


11 


08 


04:49 


3.596 


3.046 


14.6 


21.63 


± 


0.15 




2005 


Apr 


11 


08 


08:30 


3.596 


3.046 


14.6 


21.70 


± 


0.15 


176P/LINEAR 


2007 


Jan 


01 


00 


49:12 


3.162 


2.541 


16.2 


20.58 


± 


0.11 




2007 


Jan 


01 


05 


22:13 


3.163 


2.539 


16.1 


20.21 


± 


0.11 



a Heliocentric distance in AU 

b Distance from Spitzer in AU 

c Solar phase angle (Sun-object-Spitzer) in degrees 

d Expected i?-band magnitude as calculated from rotational-phase information 
inferred from the infrared data and observationally-determined H, G phase functions 
(Fig. [2} . Listed errors account for uncertainties in both rotational phase and phase 
function solutions. 



Table 2. Albedos and Radii Computed from Optical and Infrared Observations 





V 


= 0.8 a 




V 


= 1.0* 




V 


= 1.2* 




Object 


F24ixm 


r e 


PR 


F24ixm 


r e 


PR 






PR 


133P/Elst-Pizarro 


6.4 ±0.1 


1.78 


0.054 


6.5 ±0.1 


1.94 


0.045 


6.5 ±0.1 


2.09 


0.039 




6.0 ±0.1 


1.72 


0.054 


6.0 ±0.1 


1.88 


0.045 


6.1 ±0.1 


2.03 


0.039 




5.7 ±0.1 


1.67 


0.054 


5.7 ±0.1 


1.83 


0.045 


5.7 ±0.1 


1.97 


0.039 


176P/LINEAR 


10.4 ±0.2 


1.72 


0.068 


10.5 ±0.2 


1.87 


0.058 


10.5 ±0.2 


2.01 


0.050 




14.6 ±0.2 


2.04 


0.068 


14.7 ±0.2 


2.22 


0.058 


14.7 ±0.2 


2.38 


0.050 



a Assumed rj value used to compute aperture- and color-corrected 24 fj,m flux {F2n fl m) in mJy (where 
uncertainties are estimated from sky background statistics), effective radius (r e ) in km, and geometric 
i?-band albedo (pr) 
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1 "i 

i 

1 




(a) 133P/Elst-Pizarro 


(b) 176P/LINEAR 



Fig. 1. — Composite (PBCD) 24 /j,m images of (a) 133P/Elst-Pizarro (166 s total exposure time) and (b) 
176P/LINEAR (48 s total exposure time), indicated by arrows, obtained using MIPS on Spitzer. Both 
objects are point sources with no indication of cometary activity. Each panel is ~ 7.'5 by 8.'2 in size. 
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Fig. 2 — Phase function solutions for (a) 133P/Elst-Pizarro and (b) 176P/LINEAR, with best-fit IAU 
phase laws plotted as solid lines and best-fit linear phase functions plotted as dashed lines. Dotted lines 
indicate the expected range of brightness deviations from the IAU phase law due to rotation of the body. 
Observed reduced i?-band magnitudes are plotted as solid circles. 
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(e) 
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Fig. 3. — Histograms showing albedo distributions for (a) C-class (C-, B-, G-, or F-type) asteroids, where 
the superimposed black-shaded histogram only includes objects explicitly classified as C-type asteroids, (b) 
dynamical members of the Themis asteroid family, where the superimposed black-shaded histogram only 
includes those Themis members with measured albedos that have been classified as C-class asteroids (not 
all have been assigned taxonomic classes to date, however), (c) B-type asteroids, (d) active comet nuclei , 



and (e) D-type aste roids. All taxonomic classificatio ns follow the Tholen sys tem (ITholen fc Barucci 



198 



9; 



Lazzaro et aLll2004h . Comet nucleus albedos are from lLamv et al\ (2004) and B rownlee et all (|2004l ). while 
all other albedo values are fro m the IRAS Mino r Plan et Survey ( Tedesco et al.ll2004h . Objects with p < 0.075, 
designated as "cometary" bv lFernandez et all (|2005l ). are to the left of the dotted lines. 



